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Abstract 

In this paper, we contribute to the literature on energy market co-movement 
by studying its dynamics in the time-frequency domain. The novelty of our 
approach lies in the application of wavelet tools to commodity market data. 
A major part of economic time series analysis is done in the time or fre- 
quency domain separately. Wavelet analysis combines these two fundamental 
approaches allowing study of the time series in the time-frequency domain. 
Using this framework, we propose a new, model-free way of estimating time- 
varying correlations. In the empirical analysis, we connect our approach to 



the dynamic conditional correlation approach of Engle (2002) on the main 
components of the energy sector. Namely, we use crude oil, gasoline, heating 
oil, and natural gas on a nearest-future basis over a period of approximately 
16 and 1/2 years beginning on November 1, 1993 and ending on July 21, 
2010. Using wavelet coherence, we uncover interesting dynamics of correla- 
tions between energy commodities in the time-frequency space. 
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1. Introduction 



Energy commodities play an important role in economic research, as they 
affect a wide range of markets and a variety of market participants operat- 
ing on these markets. Study of the dynamics and statistical properties of 
energy commodities has become an important part of financial analysis since 
commodities became an additional tool for international diversification be- 
tween stocks, bonds and currencies. Energy-related costs play a crucial role 
in the decision making of industrial companies and entrepreneurs. A proper 
understanding of the price dynamics and interconnections between commodi- 
ties has become a fundamental concern. Moreover, energy commodity prices 
tend to have different and more extreme statistical properties than prices of 
other financial assets such as stocks, bonds, exchange rates, and correspond- 
ing derivatives. Prices of energy commodities are also influenced by various 
types of investors. Apart from standard financial investors, demand for en- 
ergy commodities consists of industrial entities and power stations using the 
commodities to produce electricity. 

The commodity markets are complex systems of interacting agents with 
different term objectives. Hence, time series resulting from this process are 
formed by a combination of different components operating at different fre- 
quencies. Standard time series econometric methods usually consider the 
frequency and time components separately. In this paper we introduce the 
wavelet approach, which allows us to study the frequency components of time 
series without losing the time information. The introduction of wavelets 
helps us uncover interactions which are hard to see using any other mod- 
ern econometric method and which would otherwise stay hidden. Moreover, 
the wavelet analysis approach is model-free. This property makes it a very 
powerful tool in comparison with other methods that rely on parameters as 
well as the estimation method. We are not the first to apply wavelets in 
the analysis of commodities. Davidson et al. (1997) proposed a form of semi- 
nonparametric regression based on wavelet analysis to study commodity price 



behavior. Yousefi et al. (2005) used wavelets for prediction of oil prices. Con- 



nor and Rossiter (2005 ) were the first to estimate price correlations by basic 



scale decomposition of time series using wavelets on commodity markets. 
They use the discrete type of wavelet transformation, while in our analy- 



sis we follow the continuous transform approach. Contrary to Connor and 



Rossiter (2005), we study the dynamics of correlations in the time-frequency 



space; our results are thus even more general. Recently, Naccache (2011) 
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used wavelets for the correlation analysis of oil prices and economic activity. 
However, that study was restricted to monthly data, with its main focus on 
longer cycles. 

In recent years, applied research on the dynamics of energy commodi- 
ties has grown significantly. Literature on the topic has been presented in 
three main categories - research on the dynamics of specific commodity time 
series, research on the macroeconomic relationships of energy commodities, 
and research on co-movements between commodities. Several studies have 



examined the topic of co-movement of commodities. For example, Lanza 



et al. (2006) examined the dynamic conditional correlations in daily returns 



of WTI oil futures prices and found these to vary significantly. Grine and 



Diko (2010) used cross-correlations to successfully construct a jump-diffusion 



model for pricing of cross-commodity derivatives of energy assets. Finally, 



Ghoshray and Johnson (2010) concluded that trends in energy prices change 



frequently and are hard to predict. 

In our research, we contribute to the discussion of co-movements and 
use a novel approach which is much easier to interpret. The application of 
wavelets enables us to study the interdependence of energy time series in 
the time as well as frequency domains, providing a deeper understanding 
of possible dependencies. The main question of our analysis is whether the 
interconnection between studied commodity markets changes significantly in 
time and varies across different investment horizons. 

More specifically, we focus on dynamic correlations based on wavelet co- 
herence between crude oil, gasoline, heating oil, and natural gas. By doing so, 
we analyze the evolution of these correlations in time as well as for different 
frequencies. This approach distinguishes between different types of investors 
with different investment horizons. Moreover, cyclical components, which 
are an important part of energy commodity dynamics, are also easily distin- 
guishable and separated between various time periods. In addition to wavelet 
coherence in the time-frequency space, we propose a new way of estimating 
model-free time-varying correlations. To relate our results to the standard 
econometric literature, we connect our approach with the standard econo- 
metric approach of Engle (2002) - dynamic conditional correlations from a 
multivariate GARCH model. 

This paper is structured as follows. After a brief introduction to wavelets, 
Section 2 presents the wavelet coherence used for the estimation of local 
correlation in the time-frequency domain, as well as the phase differences used 
to characterize the relationship between the two time series. The comparison 
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model - DCC GARCH - is also briefly introduced in this part. Section 3 
introduces the data and provides results of the empirical analysis. Section 4 
concludes by stating our main result, namely, a very interesting evolution of 
correlations between energy commodities. 



2. Methodology 



2.1. Wavelet analysis 

The wavelet transform offers localized frequency decomposition, providing 
information about frequency components. As a result, wavelets have signif- 
icant advantages over basic Fourier analysis when the object under study is 



locally stationary and inhomogeneous - see Gengay et al. (2002); Percival 



and Walden 


(2000 


); 


Ramsay 


(2002) 



In our work we use continuous wavelet 
analysis tools, mainly wavelet coherence, measuring the degree of local cor- 
relation between two time series in the time-frequency domain, and wavelet 
coherence phase differences. Before introducing wavelet coherence, let us 
provide basic definitions of wavelet and wavelet transform. 

2.1.1. Wavelet 

A wavelet is a real- valued square integrable function, ip G L 2 (IR)Q defined 



as: 



1 i. 

s \ s 



(1) 



where the term 1 / denotes a normalization factor ensuring unit variance 
of the wavelet, H'j/'u.sll 2 = 1. A wavelet has two control parameters, u and s. 
The location parameter u determines the exact position of the wavelet and 
the scale parameter s defines how the wavelet is stretched or dilated. Scale 
has an inverse relation to frequency; thus lower (higher) scale means a more 
(less) compressed wavelet, which is able to detect higher (lower) frequencies 
of a time series. In addition, there are several conditions that a wavelet 
needs to satisfy. The most important is the admissibility condition, ensuring 
reconstruction of a time series from its wavelet transform. The admissibility 
condition is defined as 



*(/) 
/ 



-df < oo, 



(2) 



1 A function x(t) is called a square integrable if x(t) 2 dt < oo. 
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where is the Fourier transform of a wavelet if>(.). The admissibility 

condition implies that the wavelet does not have a zero frequency component 
and so the wavelet has zero mean, f_ ip(t)dt = 0. Further, the wavelet is 
usually normalized to have unit energy, i.e., f_ ip 2 (t)dt = 1, implying that 
the wavelet makes some excursion away from zero. 

There is a large number of different wavelets. Each wavelet has its specific 
characteristics and is used for different purposes - see |Percival and Walden 



(2000 


); 


Adisson 


(2002 



as: 



71 



1/4 



-t 2 /2 



(3) 



Parameter w denotes the central frequency of the wavelet. We set uo = 
6, which is often used in economic applications - see for example |Aguiar- 



Conraria et al. (2008); Rua and Nunes (2009). The Morlet wavelet belongs 



to the family of complex or analytic wavelets, hence this wavelet has both 
real and imaginary parts, allowing us to study both amplitude and phase. 

2.1.2. The continuous wavelet transform 

The continuous wavelet transform W x (u,s) is obtained by projecting a 
specific wavelet onto the examined time series x(t) e L 2 (IR), i.e., 



W x (u,s) 



(4) 



An important feature of the continuous wavelet transform is the ability to 
decompose and then subsequently perfectly reconstruct a time series x(t) G 
L 2 (R): 

W x (u, s)if> UjS (t)du 



x(t) 



C 



ip Jo 



ds 

S 2 



s > 0. 



(5) 



Furthermore, the continuous wavelet transform preserves the energy of the 
examined time series, 



x 



C 



ip Jo 



\W x (u,s)\ du 



ds 

72 ' 



(6) 



We use this property for the definition of wavelet coherence, which measures 
the size of the local correlation between two time series. 
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2.1.3. Wavelet coherence 

To be able to study the interaction between two time series, we need to 
introduce a bivariate framework called wavelet coherence. For the proper 
definition of the wavelet coherence, we need to introduce the cross wavelet 



transform and cross wavelet power first. Torrence and Compo (1998) defined 



the cross wavelet transform of two time series x(t) and y(t) as 



W xy (u, 



W x (u,s)W*(u, s) 



(7) 



where W x (u,s) and W y (u, s) are continuous wavelet transforms of x(t) and 
y{t), respectively, u is a position index, and s denotes the scale, while the 
symbol * denotes a complex conjugate. The cross wavelet power can easily be 
computed using the cross wavelet transform as \ W xy (u, s)\. The cross wavelet 
power reveals areas in the time-frequency space where the time series show a 
high common power, i.e., it represents the local covariance between the time 
series at each scale. 

The wavelet coherence can detect regions in the time-frequency space 
where the examined time series co-move, but do not necessarily have a high 



common power. Following the approach of Torrence and Webster (1999), we 



define the squared wavelet coherence coefficient as: 



R 2 (u,s) 



is-^w^u,*))! 5 



S(s-i\W x (u,s)\*)S(s-i\W y (u,s)\2y 



(8) 



where 5* is a smoothing operatoiQ The squared wavelet coherence coeffi- 
cient is in the range < R 2 (u,s) < 1. Values close to zero indicate weak 
correlation, while values close to one provide evidence of strong correlation. 
Hence, the squared wavelet coherence measures the local linear correlation 
between two stationary time series at each scale and is analogous to the 
squared correlation coefficient in linear regression. 

Since the theoretical distribution for the wavelet coherence is not known, 
we test the statistical significance using Monte Carlo methods. In the testing 



Compo (1998). 



procedure, we follow the approach of Grinsted et al. (2004) and Torrence and 



2 Without smoothing, the wavelet coherence equals one at all scales. Smoothing is 
achieved by convolution in both time and scale. The time convolution is performed with 
a Gaussian window, while the scale convolution is done with a rectangular window - see 



Grinsted et al. 



(20041. 
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The use of wavelets brings with it the difficulty of dealing with boundary 
conditions on a dataset with finite length. This is a common problem with 
any transformation relying on filters. In our paper, we deal with this problem 
by padding the time series with a sufficient number of zeroes. The area 
where the errors caused by discontinuities in the wavelet transform cannot 
be ignored, i.e., where edge effects become important, is called the cone of 
influence^] ( [Grinsted et aTj |2004[ ). The cone of influence lies under a cone 



which is bordered by a thin black line. 
2.1.4. Phase 

To complete our analysis, we also use wavelet coherence phase differences 
showing us details about the delays in the oscillation (cycles) between the two 
time series under study. Following Torrence and Webster| ( |1999 ) we define 



the wavelet coherence phase difference as: 

^("• s) = tan U{s(.-w,( u ,.))}J- () 

Phase is indicated by arrows on the wavelet coherence plots. A zero phase 
difference means that the examined time series move together. The arrows 
point to the right (left) when the time series are in-phase (anti-phase) or 
are positively (negatively) correlated. Arrows pointing up means that the 
first time series leads the second one by 90°, whereas arrows pointing down 
indicates that the second time series leads the first one by 90°. Usually we 
have a mixture of positions, for example, an arrow pointing up and right 
means that the time series are in phase, with the first times series leading 
the second one. 

2.2. DCC GARCH analysis 

We would like to provide a connection of this non-traditional analysis to 
standard econometric tools. In fact, we step back from the time- frequency 
domain into the time domain to provide this connection. For this purpose, we 
compute local time- varying correlations from the wavelet coherence and com- 
pare them with the well-established multivariate concept of Dynamic Condi- 
tional Correlation Generalized Autoregressive Conditional Heteroscedasticity 



3 The cone of influence is highly dependent on the type of wavelet used - see 



and Compo (1998) 



Torrence 
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(DCC GARCH) of Engle (2002). In this part, we will provide a very basic 
overview of the model. 

The DCC estimator is a logical extension of the constant conditional cor- 
relation (CCC) model introduced by Bollerslev (1990). In Bollerslev's model, 
correlation matrix R is constant: H t = D t RD t , where D t = diag{^Jhij} and 
h ijt represents the i-th univariate (G)ARCH(p, q) process, and % = 1, . . . ,n 
at time t = 1, . . . ,T. Engle (2002) allowed R to vary in time t, thus 

H t = D t R t D t . (10) 
The correlation matrix is then given by the transformation 

R t = diag(y/qii >t , ^/q nn ,t )Qtdiag( y/q~^~t, y/Qnn,t), (U) 
where Q t = (q ij)t ) is 

Qt = {l-a-P)Q + oah-ntt-i + PQt-i, (12) 

where r\ t = y/hi~t are the standardized residuals from the (G)ARCH 
model, Q = T _1 ^rjtVt is a n x n unconditional variance matrix of rjt, and 
a and (3 are non- negative scalars such that a + (3 < 1. 

To estimate DCC GARCH, we use the standard procedure proposed by 



Engle (2002). 



3. Empirical results 

3.1. Data description 

To study the dependence between energy markets, we use the main com- 
ponents of the energy sector according to the Continuous Commodity Index 
constructed by the Commodity Research Bureau. Namely, we use crude oil, 
gasoline, heating oil, and natural gas on a nearest-future basis. The energy 
prices were collected on a daily basis over a period of approximately 16 and 
1/2 years beginning on November 1, 1993 and ending on July 21, 2010. Alto- 
gether, the sample includes 3573 daily prices for each commodity. The data 
were obtained from the Pinnacle Data Corporation. Figure [T] shows the nor- 
malized plots of the prices and Table [2] in Appendix A provides descriptive 
statistics for the logarithmic returns used in the analysis. 

In Figure [T] we can directly observe that heating oil and crude oil are 
closely related. Gasoline departs from this relationship but still has some 
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Figure 1: Normalized prices of heating oil, gasoline, natural gas, and crude oil. 

common trends during several specific periods. Natural gas seems to have 
the weakest relation to the other commodities. Let us explore these depen- 
dencies in more detail. We will first use the simple unconditional correlation 
coefficient of the returns for the whole data set, and then we will look at the 
evolution of the local correlations using wavelet coherence and connect the 
wavelet coherence results using the standard DCC GARCH approach. 

Unconditional correlations in Table [T] show that heating oil, crude oil, and 
gasoline are closely related, providing high correlations. Heating oil has the 
strongest relationship with crude oil. Gasoline is related to crude oil more 
strongly than to heating oil, but all three pairs show positive correlations of 
around 0.6. The dependence of natural gas on the rest of the group is the 
weakest, from the point of view of the correlations. 

3.2. Evidence from the wavelet coherence 

Unconditional correlations provide evidence of strong dependencies in the 
energy commodities. However, the time span of our study is quite long, so 
it may be interesting to see how the correlations develop in time. It may 
also be of interest to learn whether the dependencies vary across different 
frequencies, i.e., if there are stronger dependencies in the longer or shorter 
investment horizons. The wavelet coherence approach will be used as a tool 
allowing us to study the dependence in time as well as frequency domains. 

To assess the statistical significance of the local correlations in the time- 
frequency space, we use Monte Carlo simulations. Figure [2] shows the esti- 
mated wavelet coherence and the phase difference for all examined pairs of 
indices from scale 1 (one day) up to a scale of 256 (approximately one market 
year). Time is shown on the horizontal axis, while the vertical axis refers to 
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Figure 2: Wavelet coherence of heating oil, gasoline, natural gas and crude oil pairs. The 
horizontal axis shows time, while the vertical axis shows the period in days. The warmer 
the color of a region, the higher the degree of dependence between the pair. 

frequency; the lower the frequency, the higher the scale. The wavelet co- 
herence finds the regions in time-frequency space where the two time series 
co- vary. Regions inside the black lines plotted in warmer colors represent re- 
gions with significant dependence. The colder the color is, the less dependent 
the series are. Cold regions outside the significant areas represent time and 
frequencies with no dependence in the commodities. Thus, we can clearly 
see both the frequency and the time intervals where the commodities move 
together significantly. A continuous wavelet transform at any given point 
uses the information of neighboring data points, so areas at the beginning 
and end of the time interval should be interpreted with caution, as discussed 
in previous sections. This is also the reason we include only scales up to 256. 
From analysis of the wavelet coherence, we can observe very interesting 
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results. A first glance confirms our findings from the previous analysis. The 
three commodities of heating oil, gasoline, and crude oil show the strongest 
dependence. Heating oil is strongly related to crude oil over many periods 
and at many frequencies. Gasoline is also quite strongly related to crude 
oil. The dependence of heating oil and gasoline turns out to be restricted 
to only some periods. For example, around 1998, and 2001, these two were 
very closely related at almost all frequencies, but around the year 2000, there 
seems to be no dependence at all. This interesting finding may have its origins 
in the period of recessions with falling prices. The result will be discussed in 
the following text in greater detail. 

We can thus see that the dependence is highly dynamic as it varies in 
time. Phases, represented by arrows, do not provide any additional value 
to the analysis, as they point to the right most of the time, meaning that 
the significant local correlations are positive, but no commodity is leading 
(affecting the other one). There are short periods with changing phases and 
where one commodity seems to be leading the other, but these are not consis- 
tent so we cannot conclude that there is any directional influence. When we 
compare these three indices to natural gas, we can see very small significant 
of dependence. 

3.3. Connection to the standard analysis in time domain 

The wavelet coherence measures co-movement of two time series in the 
time-frequency space. Its main advantage is the ability to decompose the 
time- varying co-movement into different investment horizons. Since the wavelet 
coherence can be interpreted, in some sense, as a measure of local correlation, 
we would like to connect the obtained results to standard econometric anal- 
ysis, usually performed in the time domain only. In this section, we connect 
the time-varying co-movement (estimated with the wavelet coherence) with 
the DCC GARCH model, which we borrow as an representative of standard 
methods^ 

To connect the two approaches, we have to reduce the wavelet coherence 
to the time dimension only. More precisely, we compute the correlation 
coefficient based on the squared wavelet coherence coefficient for each time t 
by simply averaging wavelet coherence (those not significant at a 95 % level 



We perform only a simple-estimate comparison: this paper is not meant to provide a 
detailed discussion of the differences and capabilities of these models. 
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Heating oil 



Gasoline 



Natural gas 



Gasoline 



Unconditional 

DCC 

WTC 



0.568 

0.638 (0.137) 
0.541 (0.180) 



Natural gas 



Unconditional 

DCC 

WTC 



0.147 

0.172 (0.0235) 
0.268 (0.130) 



0.083 

0.101 (9.188xl0~ 12 ) 
0.181 (0.123) 



Crude oil 



Unconditional 

DCC 

WTC 



0.643 

0.749 (0.141) 
0.708 (0.151) 



0.612 

0.645 (0.123) 
0.626 (0.131) 



0.077 

0.109 (0.025) 
0.186 (0.129) 



Table 1: Correlation matrix - comparison of unconditional correlation coefficient with 
mean value of correlations obtained from the DCC GARCH and the wavelet coherence 
(WTC), with their standard deviations in the parentheses. 



of significance are taken to be zero) at all scales for the particular time t and 
taking the square root of this average. This reduces the wavelet coherence 
analysis to the time- varying correlations which are computed without the use 
of models. 

The squared wavelet coherence coefficients are in absolute value, hence 
negative correlations cannot be seen directly. While the wavelet coherence 
is performed with a complex wavelet, we obtain phase positions of the ex- 
amined time series which might be used for distinguishing between positive 
and negative correlations. In our particular datasets, there are no negative 
correlations (anti-phase positions of time series); hence we do not need to 
compensate for the negative correlation. 

Let us see how the time-varying correlations from the wavelet coherence 
are related to the correlations from the DCC GARCH. Table [T] provides com- 
parison of both methods to unconditional correlation. More precisely, Table 
[T] compares the mean values of the correlation using different methods for the 
whole sample. For the three most correlated pairs, the estimates of correla- 
tion from the wavelet coherence are closer to the unconditional correlations 
than the estimates from the DCC GARCH. The opposite is true for the other 
three, least correlated, pairs. In the cases with very low unconditional corre- 
lation, DCC GARCH seems to converge to a constant correlation, which is 
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(a) Heating Oil - Gasoline 



(b) Heating Oil - Natural Gas 



(c) Gasoline - Natural Gas 




Figure 3: Comparison of correlations between the DCC GARCH (dashed line) and time- 
varying correlations from the wavelet coherence (black line) for heating oil, gasoline, nat- 
ural gas, and crude oil pairs. The results are smoothed using 50-day moving averages 

close to an unconditional one. Still, except for the heating oil - crude oil pair, 
all unconditional correlations are within the standard deviation of arithmetic 
mean of correlations from the wavelet coherence as well as the DCC GARCH. 
These results point to the correctness of both approaches^] 

Figure [3] provides comparison of time- varying correlations using both ap- 
proaches^} Since we use the comparison to DCC GARCH only as an illus- 
tration of how wavelet coherence is related to standard methods, we do not 
concentrate on the differences between approaches; rather, we discuss the 
findings in the following section. 

3.4- Time-frequency dynamics of the co-movement 

The time-varying dynamics of correlations is confirmed using both ap- 
proaches. But the wavelet coherence provides a rather complex insight into 
this dynamics. Turning back to Figure [2j we can see that co-movement 
changes not only in time, but is different also at various investment horizons. 
The study of the pair of heating oil and gasoline is a good example: 



5 It should be noted that forecasting with the aid of wavelet coherences is difficult in 
the present setting due to the nature of wavelets. 
6 We use 50-day moving averages for illustration. 
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In 1998, there was a large increase in the correlation of this pair, up to 
0.8, confirmed by both DCC GARCH as well as wavelet coherence methods. 
Wavelet coherence plotted in Figure [2] shows that this strong dependence was 
present for all investment horizons from the longer ones of several months, up 
to the shorter, daily ones. In the following two years, the correlation dropped 
quickly. The wavelet coherence plot again brings the very interesting insight 
that the only dependence present in this pair around the year 2000 was 
a short-term one - up to one week. After the year 2000, the dependence 
again increased rapidly to the 0.8 level during the following year. The latest 
development of dependence among the three most correlated commodities - 
heating oil, gasoline, and crude oil - is very similar. Correlations between 
them increased rapidly to the 0.8 levels at the beginning of 2009. 

For these three most correlated commodities, the periods of high coher- 
ence around the years 1998 and 2001 are closely related to periods of reces- 
sion with falling prices. Specifically, they relate to the Asian financial crisis 
in 1998-2000, the 9/11 terrorist attacks, and the resulting fear in the mar- 
kets in 2001-2002. Finally, the current financial crisis of the period 2008-2010 
has shown a similar behavior. All three commodities again have significant 
coherences through all investment horizons. 

While we can see that the dynamics of correlation is changing rapidly 
in time, we can also see that the dynamics change in frequency as well. For 
example, the gasoline - crude oil pair shows a large dependence for investment 
horizons of several months during the years 1997-2001, while the short-term 
dependence seems to be very low during this period. The most correlated 
pair, heating oil - crude oil, shows the dynamics of the dependence in both 
time and frequency as well. The longer-term (up to six months) investors 
believe that this dependence is very high, while in the short-term investment 
horizons, the dependence varies over time much more quickly. 

We can compare wavelet coherence with the findings of entropy analy- 
sis published recently by Martina et al. (2011). The authors show that for 
crude oil, the periods with economic downturns exhibit a reduced market 
complexity in terms of reduced entropy levels. We observe similar behavior, 
because these periods are characterized by falling prices, i.e., movement in 
one direction. As a consequence, the complexity is reduced, and in the mul- 
tivariate setting the co-movement is high and significant for a large number 
of investment horizons. 



14 



4. Concluding remarks 

In this paper, we contribute to the literature on co-movement on the 
energy markets by researching the interconnections between the main com- 
ponents of the energy sector in the time-frequency space. The novelty of our 
approach lies in studying the co-movement of energy markets in the time- 
frequency space for the first time and comparing the results to the standard 
econometric tools for studying relationships between markets. 

The main finding of this paper is that some energy pairs show strong 
dynamics in co-movement in time during various investment horizons. The 
results suggest that when looking at the dependence of energy markets, one 
should always keep in mind its time-varying nature and look at it for vari- 
ous investment horizons. While the strongest dependence occurs during the 
periods of sharp price drops, it seems that the periods of recession creating 
fear in the markets imply a much higher downside risk to a portfolio based 
on these commodities. This inefficiency of the energy market is muted after 
recovery from the recession. 

Wavelet coherence also uncovered long cycles (64 to 128 days) in heating 
oil - crude oil pair that were also present in the periods outside of recession, 
or more precisely, periods of stable growth. Dependence on higher frequency 
cycles changes considerably in time. Still, the three commodities, heating oil, 
gasoline and crude oil strongly co-move, thus for the manager willing to keep 
a well diversified portfolio, the trio will imply great exposure to risk. On the 
other hand, natural gas seems to be unrelated to all three commodities for 
all investment horizons as well as the studied time periods. 

In conclusion, we uncover some interesting dynamics of the co-movement 
between energy markets in time as well as various investment horizons. Our 
findings are model-free and provide the possibility of new research on financial 
risk modeling, as they show that dynamic diversification is required in order 
to preserve a higher profit. 
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Appendix A 





Heating oil 


Gasoline 


Natural gas 


Crude oil 


Mean 


0.000368314 


0.000400265 


0.000211919 


0.000412954 


St.dev 


0.029718 


0.0320496 


0.049457 


0.0277684 


Skewness 


-1.28498 


-0.284786 


0.178664 


-0.242562 


Kurtosis 


34.5392 


6.79143 


22.049 


8.01979 


Min 


-0.520032 


-0.222545 


-0.5187 


-0.200321 


Max 


0.244974 


0.230015 


0.659246 


0.181293 



Table 2: Descriptive statistics of the daily logarithmic returns for heating oil, gasoline, 
natural gas, and crude oil in the period of November 1, 1993 to July 21, 2010. 





Coefficient 


Standard Deviation 


OlDCC 


0.04435448 


0.01306025 


Pdcc 


0.90056051 


0.03265073 


Log. Likelihood 


82605.52 





Table 3: Coefficient estimates of multivariate DCC GARCH for heating oil, gasoline, 
natural gas, and crude oil in the period of November 1, 1993 to July 21, 2010. ancc and 



Pdcc are coefficients from Equation 12 



16 



References 

Adisson, P. (2002). The Ilustrated Wavelet Transform Handbook. The Insti- 
tute of Physics, London. 

Aguiar-Conraria, L., N. Azevedo, and M. J. Soares (2008). Using wavelets 
to decompose the time-frequency effects of monetary policy. Physica 
A 387(12), 2863-2878. 

Bollerslev, T. (1990). Modeling the coherence in short-run nominal mod- 
eling the coherence in short-run nominal exchange rates: A multivariate 
generalized arch model. Review of Economics and Statistics 72, 498-505. 

Connor, J. and R. Rossiter (2005). Wavelet transforms and commodity prices. 
Studies in Nonlinear Dynamics & Econometrics 9(1). 

Davidson, R., W. C. Labys, and J. Lesourd (1997). Walvelet analysis of 
commodity price behavior. Computational Economics 11(1), 103-128. 

Engle, R. (2002). Dynamic conditional correlation: A simple class of mul- 
tivariate generalized autoregressive conditional heteroskedasticity models. 
Journal of Business & Economic Statistics 20(3), 339-350. 

Gengay, R., F. Selguk, and B. Whitcher (2002). An Introduction to Wavelets 
and Other Filtering Methods in Finance and Economics. Academic Press. 

Ghoshray, A. and B. Johnson (2010). Trends in world energy prices. Energy 
Economics 32(5), 1147-1156. 

Grine, S. and P. Diko (2010). Multi-layer model of correlated energy prices. 
Journal of Computational and Applied Mathematics 233, 2590-2610. 

Grinsted, A., J. C. Moore, and S. Jevrejeva (2004). Application of the cross 
wavelet transform and wavelet coherence to geophysical time series. Non- 
linear Processes in Geophysics 11, 561-566. 

Lanza, A., M. Manera, and M. McAleer (2006). Modeling dynamic condi- 
tional correlations in wti oil forward and futures returns. Finance Research 
Letters 3, 114-132. 



17 



Martina, E., E. Rodriguez, R. Escarela- Perez, and J. Alvarez- Ramirez (2011). 
Multiscale entropy analysis of crude oil price dynamics. Energy Eco- 
nomics 33(5), 936-947. 

Naccache, T. (2011). Oil price cycles and wavelets. Energy Economics 33(2), 
338 - 352. 

Percival, D. B. and A. T. Walden (2000). Wavelet Methods for Time series 
Analysis. Cambridge University Press. 

Ramsay, J. B. (2002). Wavelets in economics and finance: Past and future. 
Studies in Nonlinear Dynamics & Econometrics 6(3). 

Rua, A. and L. C. Nunes (2009). International comovement of stock market 
returns: A wavelet analysis. Journal of Empirical Finance 16(4), 632 - 
639. 

Torrence, C. and G. P. Compo (1998). A practical guide to wavelet analysis. 
Bulletin of the American Meteorological Society 79(1), 61-78. 

Torrence, C. and P. J. Webster (1999). Interdecadal changes in the enso- 
monsoon system. Journal of Climate 12(8), 2679-2690. 

Yousefi, S., I. Weinreich, and D. Reinarz (2005). Wavelet-based prediction 
of oil prices. Chaos, Solitions and Fractals 25, 265-275. 



18 



